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Abstract 

New aspects of a relation between lattice and dislocation structures are exam- 
ined within a physically transparent theoretical scheme. Predicted features 
originating from the lattice discreteness include: (i) multiple core disloca- 
tion structures and (ii) their dependence on the position of the dislocation 
axis. These effects, which in principle can be observed directly and may also 
manifest themselves in dislocation motion or/and transformation (cross-slip) 
characteristics, are very general and present in any crystal in which they may 
be more or less pronounced depending on the material. 
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It is widely accepted that such defects as dislocations significantly influence a number 
of properties in real materials. Thus, understanding the relationship between lattice and 
dislocation structures is one of the fundamental problems of materials physics. Despite recent 
developments of powerful atomistic simulation techniques, up to now our understanding of 
the relation between lattice and dislocation structures is based on results obtained within 
the framework of the Peierls-Nabarro (PN) model [|l]||. This model has provided both 
language for interpretation experimental/theoretical results and simple relations between 
dislocation properties and lattice discreteness characteristics (periodicity, symmetry, etc.) 

This remarkable breakthrough in understanding how dislocation properties are related 
with lattice characteristics became possible due to two features of this model: (i) its high 
tractability and (ii) its combined different length scale descriptions. 

The combined "continuum/atomistic" descriptions in the PN model follow clearly from 
the structure of the energy functionals ||, E to t, of the dislocation displacement distribution 
u(x) (here x is a distance from the dislocation axis in the slip plane) 

E tot (u(x)) = E el (u{x)) + E mi3 (u(x)) (1) 

with a linear elastic (E e i ) and a non- linear atomistic misfit energy term (E mis ) 

E mis = hY,®(u(nh-l)), (2) 

n 

where <&(u(x)) is a periodic energy profile which is often approximated by the so-called 
generalized stacking fault energy (GSF) or 7-surface ||. Indeed, the E mis term represents the 
most apparent and important lattice properties - discreteness/periodicity/symmetry - and 
allows one to investigate within the PN model the relation between lattice and dislocation 
properties. E mis can be expanded in a Fourier series as 

s=oc 2-ksI 

E mis = E° mis + h^2 J s cos — — (3) 

s=l n 

00 

where E^ is = J $(u(x))dx is independent of the position of the dislocation axis I and terms 

—00 

which are oscillatory with h being a repeat distance normal to the dislocation line and n an 
integer number that counts atomic rows in the same direction. 
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Now, a minimization of E tot {u(x)) allows one to find the equilibrium dislocation structure, 
u(x). In order to perform this minimization in analytic form, a critical approximation E mis = 
E^ is has been made JIJ. This "continuum" approximation in representing the misfit energy - 
which is supposed to describe lattice discreteness - results in an obvious inconsistency which 
has been a subject of debate for years, mainly in the context of Peierls stress determinations 
H|7|||. Thus, an interesting and fundamental issue arises - if consistently represented in the 
atomistic interaction energy, how will lattice discreteness be manifested in the structure of 
dislocations ? 

Several recent attempts to overcome this inconsistency resulted in purely numerical pro- 
cedures Thus, one of the most advantageous features of the PN model - high tractabil- 
ity and transparency - has been sacrificed. These authors also focused on the Peierls stress 
determination and demonstrated that indeed the discrete representation of the misfit energy 
brings theoretical estimates much closer to experimental results. 

In this Letter, using a physically transparent solution of the PN model with a consistent 
discrete representation for the misfit energy, we examine how lattice discreteness may influ- 
ence dislocation structure. This allows us to predict new generic features of the dislocation 
structure that are independent of the PN model assumptions, and driven by lattice discrete- 
ness such as multiple core structures and their dependence on the position of the dislocation 
axis. 

To determine dislocation structure, we perform a minimization of the total energy func- 
tional, Eq. U, with a discrete representation of the misfit energy, EqfJ, using trial functions, 
u(x), defined from the Laurent expansion PJT0|] of their derivatives pp{x) 



An J -A N Pk A 13 

ux k=m=i [x — z k ) 

where N is the maximal number, p k is the maximal order of the poles z% and AP nk are 
expansion coefficients. It is important to note that, by definition, these trial functions 
provide a minimum of E tot for an arbitrary potential (not only sinusoidal, as trial 

functions used in || to parameterize total energy functional in the convenient form) in case 



of the "continuum" approximation, E mis = E^ is @,[T(J. This choice of trial functions not 



only provides good accuracy and stability of the minimization procedure |TD] but also allows 
one to express E tot through parameters describing the dislocation structure. Indeed, the 
poles = Xj. + iu>k have a clear meaningixf gives the position and tuf gives the width 
of the partials for the screw (/3=1) and edge {13=2) components of the displacement in 
the partial cores. For example, for the ordinary dislocations dissociated into two Shockley 
partials, xf = / ± dP j2, where d is the partials separation and I gives the position of the 
whole ordinary dislocation center. In the general case with these trial functions, E tot can be 
presented as a numerical function of geometrical parameters, in particular, for an ordinary 
dislocation as a function of the set of parameters ({g} = {d, u, I}) describing the dislocation's 
structure (d,u>) and its position in the lattice (/). 

As examples, we consider ordinary dislocations for fee metal, Ir, and an ordered alloy, 
CuAu, with Ll structure. In these materials, this type of dislocation normally splits into 
two Shockley partials P,|i~0fl and represents a very typical example of dislocation structures. 
To illustrate graphically minimization of E tot , let us introduce elastic (F e i(d)) and misfit 
(F mis (d)) generalized forces which are defined as, F el (d) = —dE e i(d)/dd (the sign is chosen 
for convenience) and F mis (d) = dE mis (d) / dd (here and further we drop the (3 index since 
d l ~ d 2 ). In this definition, for a given d other geometrical parameters from the complete set 
{g} are taken to be such that they minimize E to t- Obviously, in this case the intersection of 
F e i(d) and F mis (d) gives a partial separation d which corresponds to the minimum of Eq.|l], 
provided that the second derivatives are positive. 

The generalized forces calculated according to this definition using ab-initio 7-surfaces 
(see for details) in the case of the screw orientation of the unit dislocation for Ir and 
CuAu are presented in Fig.|l] (a,c). For comparison, we also determine generalized forces 
for dislocations with simple model density distribution displacements composed of two delta 
functions, p(x) = bx8(x + d/2) + b 2 S(x — d/2) (see Fig. [I] (b, d)). In this case, we have a 
step function shaped dislocation for which E to t has a very simple form, 
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Etat = H ■ In(-) + d (5) 

and corresponding generalized forces F^ ep = H/d and = % s f, where % s f is the intrinsic 
stacking fault energy and H is a so-called prelogarithmic factor (see for example |J). Inter- 
estingly enough, for this model type of dislocation the well-known simple relation between 
the equilibrium partials separation and stacking fault energy, d = H/j is f f2[, can be easily 
recovered from the functional dependence in Eq. [|. As can be seen in Fig.[l], in the limit 
of large separation distances (d » ou), the PN model generalized forces defined within the 
"continuum" approximation approach those for the step function shaped dislocation. 

We now focus on how the oscillatory part of the misfit energy, usually neglected in the 
PN model analysis, affects dislocation structure. Remarkably enough, there are not only 
oscillations with I which can be expected from the Fourier expansion of the misfit energy, 
Eq. |3], but also oscillations with partials separation d for fixed I. As can be seen in Fig. [I], 
for Ir, despite the rather small amplitude of these oscillations, the effect of the misfit energy 
discrete representation is quite visible since the intersection of the generalized forces happens 
to be in the area where d dependence of the F e i is rather weak and so solutions are affected 
most by the F m i S oscillations. Finally, for CuAu the effect is dramatic since the amplitude 
of the oscillations is very large. It is significant that the discrete representation of the misfit 
energy not only changes quantitatively parameters of the dislocation structure (as in the case 
of Ir) but may also result in qualitatively new effects (as in the case of CuAu). As is evident 
from Fig.|], one such general effect, which is independent of the PN model approximations, 
is the appearance of multiple stable dislocation core configurations. Indeed, conclusions we 
draw in this study about the possibility of multiple core configurations are based on generic 
features of the F e i(d) and F mis (d) dependencies (see Fig.[L]) which follow from the general 



physics of the linear (elastic part) and non-linear (misfit part) lattice response |TT 



Moreover, within the proposed scheme it is possible to derive ]E2| the following convenient 
and physically transparent form for the energy functional, Eq.[l|, 

27rZ nd 

E tot = E° ot (d, u) + A(d, oj) cos — cos — (6) 



Here the / independent first term E^ ot = E e i + E^ is is the energy in the "continuum" 
approximation and the second term has an explicit oscillatory dependence both on I and d. 

This form reveals that the appearence of the above features associated with lattice dis- 
creteness are dependent on the relative contribution of the energies represented by the E^ ot 
and the oscillatory terms. In turn, the character of E® ot (d) dependence and corresponding 
generalized forces is driven by the competition of the partials attraction described by E < ^ nis 
( term which is dependent upon the 7-surface energetic characteristics and for large d, it 
can be well approximated by Eq. ||]) and the elastic repulsion {E e \ which is dependent on 
the elastic constants and in the limit d » 00, has a simple dependence, see Fig.p. The 
influence of the oscillatory term is predetermined by its amplitude A(d, u) which according 
to our analysis is strongly dependent on characteristics of the 7-surface fl3|| . 

It is important that Eq. |6| describes dislocation energetics for a wide range of d and 
uj HU and correspondingly a complex interdependence of all geometrical parameters (d, uo 
and I). Features of dislocation structure originating in this interdependence of geometrical 
parameters and their impact on dislocation energetics can be seen from the calculated d 
and I dependencies of the E tot in FigfJ. Indeed, proof that there are can be more than one 
stable dislocation core configuration (CuAu) can be seen in Fig.g(a). Next, the dependence 
of the partials separation on the position of the dislocation axis in the lattice is clearly seen 
in Fig.||](b). A comparison of the d(l) dependencies determined within the "continuum" 
approximation and with a discrete representation for E mis makes it evident that lattice 
discreteness is the origin of the variation in equilibrium dislocation structure depending on 
the position of the dislocation axis (core "relaxation"). This variation may result in changes 
in the number of stable core configurations and abrupt transitions between them (as for 
CuAu, see Fig.g(a,b)). 

These predicted features may have a profound impact on dislocation energetics. As 
can be seen in Fig.|2|(c), the dependence of the dislocation structure on the position of the 
dislocation axis may not only lower significantly the Peierls barrier (as also have been found 
in I0H) but may even modify the shape of the Peierls potential. It is remarkable that core 
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"relaxation" along with the existence of the multiple core configurations (the case of CuAu) 
adds a new feature to the Peierls potential - an additional minimum - which according to the 
model analysis may result in characteristic changes of the temperature dependence of the 
yield stress [jT6| . As can be seen clearly in FigfJ, the unusual shape of the Peierls potential 
in CuAu originates in the abrupt transitions between two stable core configurations, "1" 
and "2". 

We find that in addition to the known dislocation structure features, lattice discreteness 
is the origin of (i) multiple core configurations and (ii) their dependence on the position of 
the dislocation axis. Combination of these effects may result in rather complex variations 
of the dislocation structure over the crystal including changes in the number of stable core 
configurations and transitions between them. As a result, one may have to consider a 
distribution of the core configurations in a crystal under ambient conditions rather than one 
characteristic core structure which determines dislocation motion or cross-slip properties. 

As follows from our analysis, predicted features of dislocation structure originating in 
lattice discreteness are always present in crystals, but as we demonstrate with Ir and CuAu 
as examples, they appear more or less pronounced depending on characteristics of the given 
material. These fundamental characteristics can be identified within the proposed theoret- 
ical analysis primarily due to its tractability and physical transparency [13|]. While these 
features, namely multiple core configurations, can be directly verified, in principle, in high 
resolution electron microscopy experiments, they may also reveal themselves indirectly in 
low temperature internal friction experiments |17J and in mechanical properties which de- 



pend on elementary processes that are sensitive to the dislocation structure. Among such 
processes, we would emphasize cross-slip, where dislocation core structure and its changes 
under local stress may play an important role. 
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FIGURES 

FIG. 1. Dependence of the generalized elastic and misfit forces (in J/m 2 ) on partials separation 
d (in lattice constant units) for Ir and CuAu calculated within the PN model (left panel (a), 
(c)) and for simple model step function shaped ordinary dislocation (right panels (b), (d)). The 
misfit forces corresponding to the consistent discrete representation of the misfit energy in the PN 
model are presented by solid lines and those calculated within the "continuum" approximation by 
dotted-dashed lines. 

FIG. 2. Dislocation energy (in J/m) as (a) a function of the partials separation d, (b) a 
corresponding dependence of d on I (for more than one stable core configuration, solutions which 
are close in energy are numbered "1" and "2") and (c) a corresponding dependence of the total 
energy on the position of the ordinary dislocation center given by I in units of a repeat distance 
in the direction normal to the dislocation line (h), calculated for Ir (left panels) and CuAu (right 
panels); the ones determined within the "continuum" approximation (see text) are presented by 
the dotted-dashed lines. 
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Fig.1(a,b) O.N. Mryasov, Yu. N. Gornostyrev and A.J. Freeman 




Fig.1(c,d) O.N. Mryasov, Yu. N. Gornostyrev and A.J. Freeman 
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Fig. 2(a) O.N. Mryasov, Yu. N. Gornostyrev and A.J. Freeman 





Fig. 2(b) O.N. Mryasov, Yu. N. Gornostyrev and A.J. Freeman 




Fig. 2(c) O.N. Mryasov, Yu. N. Gornostyrev and A.J. Freeman 




